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For energy eigenfunctions of 1-dimensional tight binding model, the distribution of ratio of their 
nearest components, denoted by f(p), gives information for their fluctuation properties. The shape 
of f(p) is studied numerically for three versions of the ID tight binding model. It is found that when 
perturbation is strong the shape of f(p) is usually quite close to that of the Lorentzian distribution. 
In the case of weak perturbation the shape of the central part of f(p) is model-dependent while 
tails of f(p) are still close to the Lorentzian form. 
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I. Introduction 

Different versions of the one dimensional tight bind- 
ing model with constant near-neighbor-hopping have at- 
tracted lots of attention for years (e.g., see Al- 
though these versions of the model has the common fea- 
ture of Hamiltonian matrices being tridiagonal, they have 
been found showing distinct global behaviors such as 
localization and delocalization of energy eigenfunctions. 
For example, for the Anderson model with random di- 
agonal elements for the Hamiltonian matrix, localization 
can be established rigorously H . For the case of diagonal 
elements being quasi-periodic as in the ID quasi-periodic 
tight binding model, both localized and extended wave- 
functions have been found [^]|]-[7]] ■ Recently, another in- 
teresting version of the model has also been studied. It 
has been found that, if some correlation is introduced 
into the diagonal elements of the Hamiltonian matrix of 
the Anderson model, delocalization can appear [||-^2|. 
Clearly, global properties as localization and delocaliza- 
tion are not determined merely by the tridiagonal struc- 
ture of the Hamiltonian matrices. 

The main difference between the Anderson model and 
the quasi-periodic tight binding model lies in the global 
properties of the diagonal elements of their Hamiltonian 
matrices. More exactly, it lies in whether the diagonal 
elements are random or quasi-periodic (quasi-random). 
Since locally there is no difference between random and 
quasi-random numbers, one can expect that the differ- 
ence between local statistical properties of energy eigen- 
functions for the two models is not so distinct as that 
between their global behaviors. 

Statistical properties of eigenfunctions have been stud- 
ied extensively in recent years (for averaged shape see, 
e.g. , p3|-p9| , for fluctuations and correlations see, e.g., 
@^0|-^3j). It is known that the Random Matrix Model 
(RMM) (see, e.g., @,|5|) and the Band Random Matrix 
Model (BRMM) (see, e.g., HQ) can be used to de- 
scribe quantum chaotic systems and quasi- ID disordered 
models, respectively. Although one can expect that some 



features of the above mentioned ID tight binding model 
may be described by the BRMM with band width equal 
to 1, due to the constant near-neighbor- hopping, it is 
impossible for all its important features to be able to 
be described by the BRMM (let alone RM). Therefore, 
when studying statistical properties of the ID tight bind- 
ing model one must keep in mind its peculiar features. 

When studying statistical properties of energy eigen- 
functions the first quantity coming to one's mind may be 
the distribution of the components of the eigenfunc- 
tions. However, for the ID Anderson model, due to the 
exponential decay of \ci\ from its maximum, the distri- 
bution of Ci does not provide much information for local 
properties. On the other hand, the distribution of pi, for 
the so-called Riccati variable pi = Ci/c-i-i, is related to 
both local and global properties of eigenfunctions |2^] . In 
this paper we will study properties of this distribution. 

This paper has the following structure. In section II, 
some properties of pi are discussed. It is shown that 
f(p), the distribution of Pi, not only gives information 
for local statistical properties of energy eigenfunctions, 
but also is related to a global property, the Lyapunov ex- 
ponent. Tails of f(p) are expected to behave as 1/p 2 and 
in some cases the shape of f(p) is expected to be close to 
the Lorentzian distribution. In section III, the shape of 
f(p) is studied numerically for the Anderson model. It is 
shown that in the case of weak disorder most of f(p) can 
be fitted quite well by the Lorentzian distribution. In 
the case of strong disorder, central parts of /(p) deviate 
notably from the Lorentzian form and some features of 
f(p) can be explained by perturbation theory. Section IV 
is devoted to a discussion about the shape of f(p) for the 
other two versions of the ID tight binding model men- 
tioned above. Numerically, more or less similar results 
as for the Anderson model have been found for the two 
models. A brief conclusion and discussion are given in 
section V. 

II. Properties of pi 
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The Hamiltonian matrix of the ID tight binding model 
is of the form 



and 



Ha — ( ,6 



H,j-1) 



(1) 



Without the loss of generality, we take i, j = 0, 1, 2, • • • N. 
When TV goes to infinity one has a chain of infinite length. 
Corresponding to the above mentioned three versions 
of the model, the site energies are random numbers, 
quasi-periodic numbers, and partly random numbers, re- 
spectively, distributed in the region [— to/2, to/2]. For the 
first and third versions flat distributions will be chosen. 
Without the loss of generality, we will take w = 1 in 
numerical calculations. 

The stationary Schrodinger equation for eigenenergy 
E is 



vc l+ i + uc 4 _i + e t c t = Eci 



(2) 



with Ck, k = i — i + 1, being components of the cor- 
responding eigenfunction. Dividing both sides of Eq. (Q) 
by vci, one can see that it is convenient to introduce an- 
other quantity pi = Ci/ci-i, for which Eq. (|J) becomes 

p l+1 =a>i (3) 

Pi 



where 



E-u 



(4) 



Notice that the value of p\ is determined by Eq. for 
z = (c_x does not exist), i.e., p\ = ao- As a result, 
once the eigenenergy E is given, the eigenfunction can 
be obtained immediately from Eq. (||) for % = 1, • • • N — 
1. Eq. (||) for i = N (without pn+i which does not 
exist), i.e., pn — 1/ajv, is the condition for E to be 
an eigenenergy. Conversely, one can also calculate the 
eigenfunction from p^ = 1/cxn arid take p\ — a as the 
condition for E to be an eigenenergy. In either way, even 
for an arbitrary value of E one can calculate a sequence 
of pi from Eq. (||) (neglecting the last equation) . 

A relation between pi and p\ can be get from Eq. (0) , 



Aipi + Bi 



(5) 



where Ai and Bi are determined by recurrence relations 



Ai = a i ^ 1 A l 



A,. 



and 



Bi — at-iBi-i — Bi-2, 
respectively, with 

A 2 = a u Ai = l 



(6) 



(7) 



(8) 



B 2 = -l, Bx = 0. 



(9) 



(For more detailed properties of Ai and Bi, see Appendix 
A.) 

For a sequence of pi obtained from Eq. (^) for i from 
1 to N, when N goes to infinity, one can calculate the 
distribution of pi, denoted by f(p). For a finite N, only a 
histogram can be obtained. For the top of the histogram 
to be close to a smooth function, N must be large. For 
brevity, we also use f(p) to indicate the distribution for 
a finite sequence of pi. The distribution f(p) gives infor- 
mation for not only local fluctuation properties but also 
some global properties of eigenfunctions. For example, 
the Lyapunov exponent can be expressed as 

1 1 " 

7(2?)= lim —In 1 — 1= lim —7 m |pi| 

n^oo n Cn n— *oo n — ' 



f{p) In \p\dp. 



(10) 



As is known, for a localized eigenfunction, \ci\ usu- 
ally first increases (on average) exponentially and then 
decreases (on average) exponentially as i varies from 
to N. Let us define the Lyapunov exponent for a finite 
sequence of pi with i from j + 1 to k as 



1 1 i Cfc 
t rln — 

k-j Cj 



1 k 

— J2 ln i^ 



(11) 



i=j+l 



where [3 denotes the sequence of pi. Denoting the se- 
quence of pi of the increasing part of the eigenfunction 
as /3_|_ and that of the decreasing part as one can 
see from Eq. ^ that y(E,j3 + ) > and j(E,/3-) < 0. 
Clearly, f+(p), the distribution of pi for f3 + , and f-(p), 
the distribution of pi for /3_, can not be equal to each 
other even when iV goes to infinity. Therefore, when 
studying f(p) distributions, there are three cases that 
should be distinguished, i.e., f+(p), f-(p) and f a (p) = 
(f+(p) + f-(p)) for Pi in both (3 + and (3-. (In what 
follows we use f(p) to denote all the three cases.) 

Since it is difficult to study the shape of f{p) analyt- 
ically and directly, a related problem can be studied to 
gain some reasonable predictions for behaviors of f{p). 
Suppose one has a set of p\ with a distribution F\ (p\ ) . 
Making use of Eq. (||), one can change p\ to p 2 and 
obtain the distribution for p 2 , denoted by F2{p2)- Pro- 
ceeding this procedure, one can obtain the distribution 
for p n , F n (p n ). As will be shown in Appendix B, F n (p n ) 
does not approach any fixed distribution when n goes to 
infinity, but their tails behave as 1/p 2 and tails of f(p) 
should also decay as 1/p 2 . 

What is also worth mentioning here is that the 
Lorentzian distribution 
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f I \ a / 7r 

= {p . b) 2 +a 2 



(12) 



with 



b = — , a 
2 ' 



is invariant under the change of the variable from p to 
(oii — 1/p) when \di\ < 2. Since 



f L (p)ln\p\dp = 



(13) 



(see Appendix B), it is impossible for f±(p) to be equal 
to Jl(p) for localized states, while it is possible for f a (p)- 
However, when ^(E) is small, it is possible for f±(p) to 
be close to /l(p)- 

III. f(p) for the Anderson model 

In this section we study the form of f(p) for the Ander- 
son model by numerical calculations. First, in subsection 
A we discuss the methods for calculating f(p). Then, in 
subsection B we study the shape of f{p) for the case of 
weak disorder, i.e., when perturbation is strong. Finally, 
in subsection C the form of f(p) is discussed for the case 
of strong disorder. 

A. Methods for calculating f(p) 

There are three methods for calculating f(p) numeri- 
cally. For the first method, one can try to find out eigen- 
functions of a Hamiltonian matrix with N large enough 
for a good statistics for f(p) obtained from the data of 
one energy eigenfunction. Since when diagonalizing large 
matrices the so-called QR method usually does not give 
correct results for long tails of eigenfunctions with very 
small |cj|, we have used the following method to calcu- 
late eigenfunctions. First we chose an initial value E' 
and calculated pi,P2, ■ ■ -Pn by making use of Eq. (||) for 
i = 0, 1, ■ ■ ■ A — 1. Then, we changed E' to make pn 
satisfy Eq. (||) for i = N as exactly as possible, i.e., to 
reduce the value of err — \(j>n ~ 1/&n)/Pn\ as much as 
possible. As an example, when v = 4.0 and N = 5 x 10 5 , 
such an eigenfunction has been found with eigenenergy 
E i=s —4.0044. The value of err for this eigenfunction is 
about 10~ 14 . For this eigenfunction, \a\, on average, in- 
creases exponentially from i = to i = N. The f+ip) dis- 
tribution for this eigenfunction is given in Fig. la (circles). 
The line in Fig. la gives f-(p) for another eigenfunction 
with E w —4.002, for which |cj| decreases, on average, ex- 
ponentially from i = to i = N. This eigenfunction was 
obtained by calculating the sequence of pn,Pn-i, • • 'Pi 
and taking Eq. (^) for i = as the condition for eigen- 
solution. The value of eer for this eigenfunction is also 
about 10~ 14 . The shortcoming of this method is that it is 
difficult (if not impossible) to find out eigenfunctions for 



which both the increasing and decreasing (on average) 
parts of \a\ are large. 

For the second method for calculating f(p), one can 
diagonalize (many) Hamiltonian matrices with the same 
value of v but different realizations of the random diag- 
onal elements. There are two restrictions for the dimen- 
sion N of the matrices. One is that it must be large 
enough for eigenfunctions to be localized (exponentially 
decay). The other is that the QR method for diagonaliza- 
tion gives correct results for long tails of eigenfunctions. 
For eigenfunctions thus obtained with eigenenergies in a 
small energy region one can calculate a faip) distribu- 
tion. The shortcoming of this method is that when v 
is large, in order to obtain localized eigenfunctions the 
dimension N of the matrices must also be large, which 
requires powful computers and long calculation time. 

The third method is suggested by the fact that to ob- 
tain a sequence of pi from Eq. (|J) it is unnecessary for 
E to be an eigenenergy. In fact, when N is large enough, 
for any value of E one can calculate a /+ (p) distribution 
(7(1?) is positive). On the other hand, for p\ = l/p^-i+i 
one can get a sequence of p[ corresponding to another 
realization of the random numbers and calculate a cor- 
responding f-(p) distribution. Numerically, it has been 
found that when N is large enough the form of the f(p) 
distributions thus obtained is not sensitive to the value 
of E, that is, whether E being an eigenenergy or not 
does not have any special influence on the form of f(p). 
Therefore, for the purpose of studying the shape of f(p), 
one does not need to calculate the exact eigenfunctions. 

We have compared the above three methods numer- 
ically and found that they give the same results. Two 
comparisons between the second and third methods are 
given in Fig. 2a and 2b, respectively, for the case of 
v = 1.0 and v = 0.2. In Fig. 2a, the solid line represents 
faip) obtained by the third mothed for E — —1.0 (5 x 10 5 
Pi have been calculated) and circles give faip) obtained 
by the second mothed of diagonalizing 30 Hamiltonian 
matrices with N = 1000. For each result of the diago- 
nalizations 20 eigenfunctions with eigenenergies around 
— 1.0 have been taken for calculating the faip)- For the 
solid line and circles in Fig. 2b, E is about —0.1, 2000 
Hamiltonian matrices with N = 250 have been diagonal- 
ized and 10 eigenfunctions have been used for each ma- 
trix. The agreement between the two methods is quite 
clear in the figure. Since the third method is the most 
powerful and convenient, we will use it for the following 
calculations in this and the next sections. 

B. Weak disorder 

As mentioned in section II, f+ip) and f-(j>) must be 
different for the Anderson model, since "fiE, (3) for the in- 
creasing and decreasing parts have different signs. How- 
ever, in the case of weak disorder Lyapunov exponents 
jiE) are usually quite small. For example, when v = 4.0, 



3 



j(E ~ —4.0) is less than 1.0 x 10~ 3 . So it is possible for 
f+(p) to be close to f-(p). Indeed, Fig. la shows that 
they are quite similar. 

As discussed in section II, the Lorentzian distribution 
in Eq. ([l2]) is a fixed distribution for the mapping in 
Eq. (||) and when j(E) is small f±(p) may be close to 
Jl(p)- Therefore, it is natural to fit the f(p) distribu- 
tion in Fig. la by the Lorentzian distribution. Since a{ in 
this case is E/v fts — 1, one can expect that parameters 
b and a for the Lorentzian distribution may be close to 
cti/2 sa —0.5 and y/l — b 2 w 0.87, respectively. Indeed, 
the best fit for the dots in Fig. la by the Lorentzian distri- 
bution gives b ~ —0.53 and a w 0.85. The result is given 
in Fig. lb. It can be seen that the distribution /+(p) is 
indeed quite close to the Lorentzian form. 

For the other values of E when v = 4.0, it has been 
found that when \E\ is larger than 0.2 and within the 
region for eigenenergies (roughly speaking \E\ < 2d), the 
shapes of f± (j>) can also be fitted quite well by the same 
Lorentzian distribution as for E w —4. When \E\ be- 
comes less than 0.2, the top of f(p) will gradually de- 
viate from that of the Lorentzian form, while tails are 
still of the Lorentzian form as in Fig. lb. As an exam- 
ple, we present in Fig. 3 central parts of f+(p) for E = 
and its fitting curve of the Lorentzian form. For compari- 
son, we have also given f+(p) distributions for E — —4.0, 
E = —7.5 and the corresponding fitting curves. The dif- 
ference between f+(p) and f-(p) has also been studied 
for < \E\ < 0.2 and has been found larger than that 
for the case of \E\ > 0.2. When E = 0, f + (p) and f_(p) 
are almost the same. When \E\ becomes less than 0.2, 
the form of f a (p) will also gradually deviate from that of 
the Lorentzian distribution in the top region. 

When the value of v is increased, similar numerical re- 
sults have also been obtained for \E\ < 2v. Particularly, 
in the case of v going to infinity (or, equivalently, let- 
ting w be zero while keeping v a constant), the shape 
of f(p) can also be fitted quite well by the Lorentzian 
distribution except for some special values of E/v, e.g., 
or —1, For these values of E/v some Pi are zero and 
the mapping in Eq. @ can not procceed. But these 
special values are unimportant because they do not cor- 
respond to any eigensolution of the Hamiltonian matrix. 
Since when w — eigenfunctions are extended and the 
Lyapunov exponents are zero, the form of f(p) may be 
exactly the Lorentzian form when N goes to infinity. 

When v is decreased from 4, deviation of f±(p) from 
the Lorentzian distribution will also appear in the region 
of large \E\. As v becomes smaller the region of devi- 
ation will become larger. For example, when v = 2 for 
most of the values of \E\ < 2v the shape of f±(j>) can be 
fitted quite well by the Lorentzian form, but when v = 1 
deviation of f± (j>) from the Lorentzian form has been de- 
tected in the top region in most cases. Correspondingly 
when v = 1 obvious difference between f+(p) and f—(p) 
has also been found. Since f+(p) and f-{p) can compen- 



sate each other, f a (p) may be of a better Lorentzian form 
than f±(p). In deed, when v = 1 and 1.7 > \E\ > 0.5, 
the form of f a {p) has been found of a good Lorentzian 
form still in the top region. 

The above phenomena for relatively small values of v 
can be explained by some results of the Appendix C. As 
shown there, an eigenfunction can be divided into pertur- 
bative and non-perturbative parts. For the perturbative 
part, Eqs. (C31) and (C32) tells that pi behaves dif- 



ferently for the increasing and decreasing parts of the 
eigenfunction. Also as shown in Appendix C, when v be- 
comes smaller the perturbative part will become larger, 
while for a fixed v eigenfunctions with large \E\ have 
larger perturbative parts than those with relatively small 
\E\. As a result, the difference between f+(p) and f-(p) 
will enlarge as v becomes smaller or \E\ becomes larger. 
Since some dynamical effects remain in A a (j — ► j + 1), 
Pi of perturbative parts will also make f+(p) and f-(p) 
deviate from the Lorentzian distribution (cf. the next 
subsection) . 

C. Strong disorder 

In the case of strong disorder (for small values of v), 
considerable difference between f+{p) and f—(p) has been 
found. An example is given in Fig. 4 for v = 0.2 and 
E = —0.1, where the circles and the dashed line represent 
/+ (p) and /_ (p) , respectively. Here /+ (p) and /_ (p) are 
distributions of Pi and 1/pi, respectively, obtained from 
Eq. (|J). As discussed in section II and Appendix B, long 
tails of f(p) should behave as 1/p 2 . In Fig. 5a and 5b 
we show the above f+(p) and f-(p) distributions (dots) 
in logarithm scale and the corresponding fitting curves 
(solid lines) of the Lorentzian form for their tails. For 
the f+(p), due to the form of the central part, we can fit 
only one of its two tails by the Lorentzian distribution 
once. For the other tail similar fitting can also be done. 
It is quite clear from the figure that long tails of both 
f+(p) and f—(p) can be fitted well by the Lorentzian 
distribution. 

In Fig. 4 it can be seen that / + (p) changes only a little 
when p varies from —2 to 1. This can be explained by 
perturbation theory. In fact, Eq. ( p31 ) in Appendix C 
shows that pi for the perturbative part of an eigenfunc- 
tion satisfies 



1 



Pi 



where, according to Eq. ( p!8 ), 



A a (j-l^j) 



E-e 



j'-i 



+ 0(v 3 ). 



(14) 



(15) 



To the first order of approximation, pj « [E — e,-_i)/w. 
So when v is small, the contribution of perturbative parts 
tends to make f+{p) close to a platform in the region of 
p 6 [{E - 0.5)/v,(E + 0.5) /v}. Since there also exist 
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non-perturbative parts for which Eq. ( |14| ) does not hold 
and when v is not quite small higher order of approxima- 
tion should also be considered, the region of the platform 
should be narrower. This can be seen quite clearly in 
Fig. 6, in which f+(p) and f-(p) for v = 0.1 and E = 0.1 
are presented. 

IV. Two other versions of ID tight 
binding model 

A. Paired correlated Anderson model 

When some correlation is introduced into the diago- 
nal elements of the Hamiltonian matrix of the Anderson 
model, derealization can appear. For example, one can 
arrange e, into pairs with 62k = £2fe+i, where t2k for dif- 
ferent pairs are random numbers |lQ] . In some energy 
regions, localization lengths for this model can be far 
larger than that for the corresponding standard Ander- 
son model (for some special energies localization lengths 
may be infinite). But for local statistical properties as 
the shape of f(p), it is reasonable to expect that the 
two models may show similar behaviors. In deed, similar 
properties as shown in section III for f(p) of the Ander- 
son model has also been found numerically for f(p) of this 
model (Fig. 7). Particularly, in the case of weak disorder, 
central parts of f(p) for \E\ around zero can be fitted by 
the Lorentzian form better than for the corresponding 
case of the standard Anderson model (see Fig. 7b com- 
pared with Fig. 3c). 

B. ID quasi-periodic tight binding model 

For this model, we take £j in Eq. (|l|) as 

d = cos(2iria)/2 (16) 

where a is an irrational number. Numerically it has been 
found that when v is large, e.g., v > 4, the form of f(p) 
can by fitted quite well by the Lorentzian form (Fig. 8a) 
when E ^ and /+ (p) is almost the same as /_ (p) since 
eigenfunctions are extended. Even for quite small values 
of \E\, the fitting is also quite good. But when E = 0, 
f(p) has a high peak in the neighbourhood of p — 0. Sim- 
ilar to the Anderson model, when v is decreased from 
4, the form of f(p) will deviate gradually from that of 
the Lorentzian distribution, especially for eigenfunctions 
with relatively large \E\. This is also due to the fact 
that for these eigenfunctions the perturbative parts be- 
come larger as v becomes smaller. In some cases it has 
been found that although f+{p) and f—(p) deviate from 
the Lorentzian form notably they are quite close to each 
other. When v is smaller than 1, many localized states 
can be found. For most of these localized states f+(p) 
have shapes much more irregular than that for the An- 
derson model. In Fig. 8b we give an example for v = 0.2 
and E = —0.1. Circles represent f+(p) and the solid 



line is the fitting curve of the Lorentzian form for the far 
right tail. The far left tail and the middle right tail (for 
p between 3 and 7) can also be fitted by the Lorentzian 
form but with different fitting parameters. 

V. Conclusions and discussions 

The distribution f(p) of ratio of nearest components 
of energy eigenfunctions, pi — Ci/ci-%, is studied numeri- 
cally in this paper for three versions of the 1-dimensional 
tight binding model: the Anderson model, a quasi- 
periodic tight binding model and a paired correlated An- 
derson model. It is shown that f(p) for increasing and 
decreasing parts of eigenfunctions (denoted as f+(p) and 
f-{p), respectively) should be treated separately. Nu- 
merically the following results have been found for all 
the three models. (1) In the case of strong perturba- 
tion (weak disorder), when \E\ is not close to zero, f+(p) 
and /_ (p) are close to each other and can be fitted quite 
well by the Lorentzian distribution. (2) In the case of 
weak perturbation (strong disorder), f+(p) deviates no- 
tably from /_ (p) and only their tails can be fitted well by 
the Lorentzian distribution. Some features of the shape 
of f+(p) in the case of strong disorder for the Anderson 
model can be explained by perturbation theory. 

Eq. (||) shows that once an energy E is given, 
Pi , P2 , • • • Pn can be calculated readily from the i = 
0, 1, • • • N — 1 cases of the equation, despite of whether 
E is an eigenenergy or not (Eq. (|^) for i = N is the 
condition for E to be an eigenenergy). Numerically it 
has been found that E being an eigenenergy or not has 
no influence on whether f(p) can be fitted well by the 
Lorentzian form. That is to say, the closeness between 
the distribution f(p) and the Lorentzian distribution in 
the case of weak disorder is not specific for eigensolutions. 
It is in fact determined by the form of the mapping in 
Eq. (||). We would like to point out that this property 
of f(p) is, in fact, related to a more general problem, 
that is, some information for properties of eigenfunctions 
may be obtained without diagonalizing the Hamiltonian 
matrix. 
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APPENDIX A: PROPERTIES OF Aj AND Bi 

In order to express Ai and Bi explicitly by making use 
of aj., • • • , cti, we first introduce a quantity Sji(0 < j < i) 
defined by 



Sji = CijOlj+i 



(Al) 



Taking out k pairs of neighbouring ai arbitrarily from 
Sji, one has a product of (i — j + 1 — 2k) ay , denoted as 
S^ . Let us define another quantity Dji as 

^* = Si*-£ 5 S w) +£$ m) -"- 



(fc,m) 

Then, making use of Eq. (||) one can proof that 
Ai = Du 

B t = -D 2l 



(A2) 



(A3) 



The following relations can be readily obtained from 
Eqs. (§), (0) and @, 



D j(i-i)Pj - D (]+l)(i-l) 



Pi = 



D (i+i)i -PiD 



(j+i)(*-i) 



Djj - piDj(i_i) 



(A4) 



(A5) 



Eq. (g) is a special case of Eq. ( A4 ) . Conversely, p\ can 
also be expressed as a function of pi , 



Pi 



Ai -piAi-i 
Making use of Eq. (||) and then ( |A6| ) , one has 
P1P2 ■ ■ -Pn-i = An-xpx + B n -1 

_ A n B n -i — A n -iB n 

A n PnA n — i 



(A6) 



(A7) 



On the other hand, from Eq. (A5) for i = n and 
j = f, 2, •••,?! — f, the following relation can be obtained 
readily, 



P1P2 ■ ■ -Pn-i 



A n PnA n —l 



(A8) 



Then, substituting Eq. (A8) into (A7) one has 

A n B n-l - An-lB 7l = 1. (A9) 



So for large \A n \ and \B n \, A n /B n is approximately a 
constant. 

An interesting result of Eqs. ( |A9| ) and (@) for the case 
of large \A n \ is that different p\ gives approximately the 
same p n . This means if \A n \ approaches infinity when n 
goes to infinity, the distributions f(p) for sequences start 
from different p\ must be the same. 

In principle, once Ai and Bi are known, explicit expres- 
sions for properties of the corresponding cigcnfunction 
can be obtained. For example, from Eqs. ([l(]) and (||), it 
can be seen that the Lyapunov exponent is determined 
by properties of A n and £?„, 



7 (£) = lim -\n\A n + ^\ 

n— >oo n pi 



(A10) 

However, to obtain practically useful expressions for Ai 



and Bi from Eqs. ( A2) and (A3) is mathematically quite 
difficult. 



APPENDIX B: F N (P) AND F L (P) 

For a set of p\ with some distribution F\(j>\), from 
probability theory (see, e.g., p9j ) and Eq. (||) it is easy 
to find out the relation between F n (p n ) and Fi(pi) 



F n (p n ) = [pxPi ■ ■ ■ Pn-if Fi(pi)\ Pn 



(BI) 



where p\,p2,- • ■ ,Pn-i on the right hand side are func- 
tions of p n - From Eqs. (A6) and (A8), it can be found 
out that 



pB n 



B rl 



( A n - pA n - 1 ) 2 An- pA n _ i 



(B2) 



where for simplicity the subscript of p n has been omitted. 
Due to the fact that ai is not a constant, F n (p) can not 
approach a fixed distribution as n goes to infinity. 
The sum of F n (p) , 



1 n 

F(p,n) = ~J2 F i(p)> 



(B3) 



is more important than F n (p), since it is just the average 
of the distributions f(p) for pi obtained from the set of 
Pi by Eq. (||). Numerically it has been found that when 
f(p) is close to /l(p) in Eq. ([[2]), F(p, n) first approaches 
fhip) when n increases from n = 1, despite of which 
form Fi(pi) is. However, when n becomes larger than 
some quantities F(p, n) will show quite large fluctuations 
if the Lyapunov exponent is positive. This is not difficult 
to understand since positive Lyapunov exponent makes 



\A n \ increases exponentially and Eq. (B2) shows that 
F n (p) has peaks in a neighborhood of q n = A n /A n -i 
which becomes smaller as \A n \ increases. Equation (B2) 



also shows that the form of F\ has considerable influence 
on the form of F n (p) only in the neighborhood of q n . 



G 



Out of the neighborhood, F n (p) behaves like l/(p — q n ) 2 . 
Therefore, F(p, n) (so /(p)) should behave as 1/p 2 in the 
tail region. 

In order to calculate the value of 7 — J fi,(p)ln\p\dp 
for fh{p) in Eq. (|l2|), we first change the variable to 
q — on — 1 1 p., which gives 



a/ir 



l-oc (l - b ) 2 + a 2 a i - 1 
Then, changing q to x — 2b — q we have 

a/ir 



\dq. 
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(a; + b) 2 + a' 



: ln\x\dx. 



(B4) 



(B5) 



Finally, after changing x to y = —x one can see that 
7 = —7 which gives 7 = 0. 



APPENDIX C: PERTURB ATI VE STUDY FOR ID 
ANDERSON MODEL 

In this appendix we study ID Anderson model by 
making use of a generalization of the so-called Brillouin- 
Wigner perturbation expansion (GBWPE) H|. Accord- 
ing to GBWPE, when perturbation is not extremely 
strong, each eigenfunction can be divided into two parts, 
one of which can be expressed as a convergent pertur- 
bation expansion by making use of the other. For the 
former part, an expression for pj will be given in this ap- 
pendix, which is useful for understanding some features 
of the shape of f(p) when perturbation is not strong. 

For the ID tight-binding model with Hamiltonian H = 
H° + V, we denote eigenstates of H° as \k > and eigen- 
statcs of H as la >, 



H°\k >- 



El\k > 



H\a >= E a \a > 



(CI) 



(C2) 



(E% = £fe in Eq. (|])) For the sake of completeness, we 
first give a brief discussion for the above mentioned GB- 
WPE. Dividing the set of basis states \k > into two parts 
denoted by Sp (including Np basis states) and Sq (in- 
cluding Nq basis states), respectively, one has two pro- 
jection operators 



P= £ |fc><fc|, 



= 1 — P. 



(C3) 



\k>£S p 



Subspaces related to P and Q will be called in the fol- 
lowing the P and Q subspaces , respectively. Split- 
ting an arbitrary eigenstate |a > into two orthogonal 
parts |a >= \t > +\h > where \t >= P\a > and 
\h >= Q\a >, it can be shown, by making use of the 
stationary Schrodinger equation, that 



a >= \t > + 



E n - H° 



QV\a > 



The iterative expansion of Eq. ( C4 ) gives 



a >= t > + 



QV\t>+{ 



1 



E a -H° 



(C4) 



QV) 2 \t > 



+ ( 



E n -H a 



QVf\t> 



if 



where 



lim <T£\TZ>=0 

71 — ^OO 



\TS >^ ( 



1 



E a -H° 



QV) n \a > 



(C5) 



(C6) 



(C7) 



Here the eigenvalue E a has been treated as a constant. 



Eq. (C5) is just the above mentioned generalization of 
Brillouin-Wigner perturbation expansion (GBWPE). For 
convenience, in what follows let us use the following no- 
tations: \i > denoting a basis state in the P subspace 
and I j > denoting a basis state in the Q subspace. Then, 
Eq. p3) gives 



-;<j\a >=X]( 



lEP 



E n 



E 



v jk 



v kl 



k,l£Q 



v. ik 



v kl 



E°E a 



E° k E a 



E" 



■)C ai (C8) 



Generally to say, < T™\T™ > vanishes as n — > 00 if 
P and Q subspaces are such chosen that, for any state 
\j > in the Q subspace, \E a — E®\ is large enough com- 
pared with V, i.e., if all the basis states \k > with small 
\E a ~ E® \ are in the P subspace. So there are generally 
many P (and co rres pondingly Q) subspaces ensuring the 
validity of Eq. (C5). In these P subspaces there exists 
one with the minimum number of basis states. Projection 
operators related to this P subspace and the correspond- 
ing Q subspace are denoted by P m and Q m , respectively. 
\t a >= P m \ a > wm be called in the following the non- 
perturbative part of |a >, and \h a >= Q m |a > the per- 
turbative part, which can be expressed in terms of \t a >, 
E a , V and H° as shown in Eq. (C5). In what follows, 
we will discuss the P m and Q m subspaces only and for 
the sake of brevity we will omit the subscript m. It is 
clear from the above discussion that the subspace P will 
increase (correspondingly, Q will decrease) as the pertur- 
bation v becomes stronger. It can also be seen that for 
the Anderson model eigenfunctions with larger \E — eTj 
have larger perturbative parts. 

In order to study C a j of the perturbative part of |a > 
in Eq. (C8), we make use of the concept of path, in 
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analogy to that in the Feynman's path integral theory 
p0| |. For (n + 1) basis states |fco >, \ki >,---|fc„_i > 
, \k n > with \ko >,---|fc„_i > in the Q subspace and 
\k n > in either the Q or the P subspace, we term the 
sequence fco —> k± — ► • • • — > fc n _i — > fe„ a pai/i of n 
paces from feg to k n , if Vjtfe' corresponding to each pace 
is non-zero. Clearly paths from ko to k n are determined 
by the structure of the Hamiltonian matrix in the H° 
representation. Attributing a factor Vkk'/(E a — E%) to 
each pace k — > k' , we define the contribution of a path 
s (from j to i) to C a j, denoted by f as (j — > i), as the 
product of the factors of all its paces. Then, C a j in Eq. 
(|C8|) can be rewritten as 



ieP 



where 



A a (j -> i) = ^ fasti 



(C9) 



(CIO) 



|a > is just this basis state \k >. Let us denote it by 
\i >. Then, for each \j > in the Q subspace, according 
toEq. @ 



C a 



A a (j * i)Cai 



(C15) 



Due to the specific form of perturbation as shown in Eq. 
([j]), paths starting from j < i can never reach points 
k > i and vice versa, i.e., paths can not cross "z". From 
the discussion of Eqs. ( Cll ) and ( C14 ) , we know that 
C a i given by Eq. ( |C15| ) satisfy Eq. ( Cll ) naturally and 



the value of E a is determined by the equation 

vCd-x + E°C m + vC ai+1 = E a C m (C16) 



where C a i-i and C a i+i, given by paths left to and righ t 
to j, resp ectively, are functions of E a (see Eqs. ( C15 ) 
and (|C1C|)). 

For simplicity, let us first discuss the case j < i. Ac- 
cording to the definitions of A a (j — > i) and fasti ~ * *)> 
one can find out that 



and s denotes possible paths from j to i. 

An interesting property of C a j given by Eq. (C£) is 
that they satisfy equations 



E H JkCc 



k E a Caj 



(Cll) 



Aa(j - *) = A «ti - J + + !-»*), (C17) 

where 

-» J + 1) = /a., (j -» j + 1) + /as 2 (j -» j + 1) 



even for arbitrary values of E a a nd C a j. In fact, from 
the definition of A a (j ~ > *) m Eq. ( |C10| ), it can be shown 
that 



A a (j 



L .E ol -Ef 



E a - Ei 



Substituting Eq. QC12| ) into Eq. (JC9D, one has 



E 



^ ~^.E a - E* 



aj' 



E 



E a -E° 



(C12) 



(C13) 



which gives Eq. (Cll). Therefore, the values of E a and 



C a i must be determined by the other part of the eigen- 
equation (p2j), i.e., by equations for C a i, 



(C14) 



and the normalization condition (except for a common 
phase) . 

Now let us apply the above results to the study of the 
Hamiltonian in Eq. (|l|). Firstly, we discuss the case in 
which the value of v is so small that there is only one 
basis state \k > satisfying \E a — E®\ < v/2. (For such a 
small perturbation strength Brillouin-Wigner perturba- 
tion expansion is also convergent.) In this case, the P 
subspace corresponding to the non-perturbative part of 



+ fas 3 (j - 3 + 1) + UsA3 -> 3 + 1) + 0(v 5 ). (C18) 



Paths si, S2, S3 and S4 in Eq. flClq) are si : j — ► j + 1, 
s 2 : j -> j + 1 -> j -> j + s 3 : j -> j + 1 -> j + 2 -> j + 1, 
and s 4 : j '• ^ j — 1 — » j — » j + 1, respectively, and 0(d 5 ) 
represents the contribution of paths with more than 3 
paces. Substituting Eq. (C17) into (C15), one has 



a 



pj+i 



1 



(C19) 



and 



C Q j = Aa(j ~> J + ^)Aa{] + 1 -» J + 2) 



Thus, on average, 



(C20) 



(C21) 



wher e A i s the average value of A a (j ~ * 3 + !)• From 
Eq. (|C18| ) one can see that |A| < 1, so |C aj -| increases 
exponentially as j approaches i. When j > i, similar 
arguments lead to the following relations 

Aa(j -> *) = ^a(j -» J - l)^a(j ~ 1 ~> »)» ( C22 ) 



Caj — 1 



A a (j^i-1) 



(C23) 
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and 



;(m) 



Caj ~ (4"* *Cq 



Secondly, we discuss the case of v being larger than 
only a few \E a — E%\. (This is a case in which Brillouin- 
Wigner perturbation expansions diverge.) In this case, 
the P subspace corresponding to the non-perturbative 
part of | a > is composed of more than one basis states 
with E® close to E a distributed randomly in the region 
[—w/2, to/2] . Since N p , the number of basis states in the 
P subspace, is much smaller than N, the total number 
of basis states, labels i are generally separated by labels 
j of | j > in the Q subspace and there is generally no 
successive labels of i. On the other hand, Sq, the set of 
\j >, is also separated into a series of subsets by \i >. For 
convenience, we denote the basis states \i > in increasing 
order by |iW >, |i( 2 ) >, . . . >, . . . |j(^) > with iW < 
j( 2 ) < . . . < . . . < i( N p) and denote the corresponding 
series of subsets of S Q by S { q\s q 2 \- ■ ■ S ( Q q) , • • •. 

To apply GBWPE to this case of perturbation 
strength, let us consider a state \j > in Sq'' satisfying 
j(<?) < j < . Due to the same reason as in the first 

case, paths starting from j can not cross i^ q ' and i^ q+1 \ 
that is, they are restricted in the region \jS q \v- q+1 }]. 
Then, according to Eq. ( |C9| ) 

C aj = A a (j -> *W) ■ C Qi( „ + A Q (j -> ■ C m<9+1 , . 

(C25) 

Similar to Eq. (|C17|) in this case we have 

4*0 = 4,0 - j - 1)4,0' - 1 - * {q) ) 

A a (j - z (9+1) ) = 4,0 -> j + 1)4,0' + 1 

(C26) 

Then, 

C QJ - = 4,0 -» j - 1) • •■ A a (i^ + 1 -> i^C^,, 

+ 4,0 - j + 1) •• • 4(* (9+1) - 1 - ^ +1) )C„(, + D 



(C27) 



(4V 1 " c ai (,) + (4* " 3 c m ( q +i) 



Since only a few \E a — E® \ for \i > in the P subspace are 
smaller than v, Eq. ( |C1§| ) also gives \A\ < 1. 

There are iVp equations as Eq. (C14) which determine 
the value of E a and give (Np — 1) relations between 
C a i. What is of special importance here is that E a is 
determined mainly by one of these equations only. In 
fact, suppose \C ai { m ) | is the largest one among \C a i\- For 



1 and j' m 



Am) 



1, Eq. flC27D shows that 



(C24) n 



both C a j m and C a j/ are mainly given by terms including 



From the discussions for Eqs. ( C14| ) and ( C16 ) 



one can see that E a is mainly determined by equation, 



E a — E 



vA a (j -> iW) + «4»(j" - l(m) )> (C28) 



i.e., by contribution of paths around i^ m >. 

As a result, for any other state \i > with i — i( q > ^ 
i^ m \ it is generally impossible for both C a j and C a y 
(j = i — 1 and j' = i + 1) to be mainly determined by 
terms including C a j. In fact, if that is possible, one will 



have an equation similar to Eq. ( C28 ) , which means that 



E a is mainly determined by the contribution of paths 
aroun d i. Generally, to say, this contradicts with Eq. 

Then, at least one of C a j and C a ji should have a con- 



siderable contribution from the other C Q i< (i 1 



£(«-!) 



,•(3+1)1 



For example, Suppose C a j has a consider- 



able contribution from C a i> with i' — l \ This means 
C a j is in the same order of magnitude as 

(Af-t'-Vc^ 



(see Eq. ( p27 )). Since according to the eigen-equation 



vC a 



(E a — E®)C a 



vC a 



(C29) 



one can see that C a i is nearly in the same order of mag- 
nitude as C a j. Thus, for % satisfying i > % > i' Eq. 



27h gives 



a 



(Ay 



Co 



(C30) 



Similar result can also be obtained if C a j' has a con- 
siderable contribution from C ai { q +\) . Therefore, when k 
changes from v q ' to i' 9+1 ) or from to v q ~ 1 ', \C a k\ 
should on average increase exponentially in at least one 
of the two cases. Specifically, if = 1 for q = 1, 
\C a k\ should on average increase exponentially when k 
increases from to . Similarly, for the case of 
j(JVp) _ TV, \C a k\ should increase exponentially when k 
decreases from N to i^p -1 ). 

Now it is easy to show that the eigenf unction of \a > 
should be localized. Let us start from the basis state |1 >. 
If 1 1 > is in the P subspace, as discussed above, \C a k\ 
should on average increase exponentially when k changes 
from 1 to i {2 \ If |1 > is in the Q subspace, Eq. flC27|) 
shows that \C a k\ also on average increases exponentially 
until k reaches the first i = . Then, making use of 
the result of the above paragraph, it can be found that 
|C a fc| should increase exponentially as k changes from 1 
to i'" 1 ). Similarly, it is easy to see that |C a fe| should on 
average decrease exponentially when k changes from i*-" 1 ' 
to N. Thus, the eigenfunction of \a > has the property 
of exponential decay. Applying this result to Eqs. (C25) 
and flC2q), one can found out the following relations 
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1 



p.! 



and 



pj w A a {j 



(C31) 



(C32) 



for the increasing and decreasing parts of eigenfunctions, 
respectively. 

Thirdly, let us increase the value of v further. In this 
case, there will be some successive labels i. That is to 
say, some of the subsets of Sp separated by Sq may have 
more than one basis states \i >. When the numbers of 
basis states \i > in the subsets of Sp are small compared 
with the numbers of basis states \j > in the subsets of Sq , 
following similar arguments as in the second step case one 
can see that when \A\ < 1 energy eigenfunctions should 
also be localized (exponential decay). Finally, when the 
number of basis states in Sp is larger than that in Sq, 
exponential decay of eigenfunctions can not be proved 
by the method used above, but can be proved by other 
methods ( f§). 

However, in these two cases the method is still use- 
ful for studying properties of Pj and two results can be 
obtained. First, for perturbative parts in a region of 
j € [Ji , J2] , if one of the values of |C a7l | and |C QJ2 | is 



far larger than the other, (C31) and (C32) can still be 
obtained from equations similar to Eqs. ( p25j ) and ( C26| ) 
Second, for \E a \ > w/2, since larger \E a \ means larger 
\E a — e,|, eigenfunctions with larger \E a \ have larger per- 



turbative parts and so have more p.,- satisfying Eqs. (C31) 
or flC3p . 
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FIG. 1. (a) A comparison (in logarithm scale) between 
the distribution ,f+(p) for an increasing eigenfunction with 
E ~ —4.0044 (circles) and the /_ (p) for a decreasing eigen- 
function with E « —4.002 (solid line) when v = 4.0 and 
N — 5 x 10 5 for the Anderson model, (b) f+(p) in (a) and 
the fitting curve of the Lorentzian form 



FIG. 2. Comparisons between f a (p) obtained by the 
second method of diagonalizing many Hamiltonian matri- 
ces (circles) and that by the third method of making use 
of Eq. (3) (solid lines) for (a) v = 1.0, E = -1.0 and (b) 
v = 0.2, E = -0.1. 
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FIG. 3. Central parts of f+(p) (circles) for the Anderson 
model and the fitting curves (solid lines) of the Lorentzian 
form for v = 4.0 and (a) E = -7.5, (b) E = -4.0, (c) E = 0. 

FIG. 4. Central parts of f+(p) (circles) and f-(p) (dashed 
line) for v = 0.2 and E = —0.1 for the Anderson model. 

FIG. 5. (a) The fitting curve (solid line) of the Lorentzian 
form for the left tail (dots) of the /+ (p) in Fig. 4 (in logarithm 
scale) . (b) Same as (a) for tails of the /_ (p) in Fig.4. 



FIG. 6. Central parts of f+(p) (circles) and f-(p) (dashed 
line) for v = 0.1 and E = 0.1 for the Anderson model. 



FIG. 7. Central parts of f+(p) (circles) for the paired 
correlated Anderson model and the fitting curves (solid lines) 
of the Lorentzian form for v = 4.0 and (a) E = —4.0, (b) 
E = 0. 



FIG. 8. f+(p) (circles) for the ID quasi-periodic 

tight binding model and the fitting curves (solid lines) 
of the Lorentzian form. (a) v = 4.0, E = —4.0, (b) 
v — 0.2, E — —0.1 and the fitting curve is for the far right tail 
only. 
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Fig.5 




